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Abstract. The higher order supersymmetric partners of a stationary periodic potential are stud- 
ied. The transformation functions associated to the band edges do not change the spectral struc- 
| ture. However, when the transformation is implemented for factorization energies inside of the 

forbidden bands, the final potential will have again the initial band structure but it can have 
bound states encrusted into the gaps, giving place to localized periodicity defects. 
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Nowadays there is a growing interest in constructing exactly and quasi-exactly solv- 
able potentials which could serve as models in various physical situations (see e.g. [T| 
and references therein). There exist simple generation techniques, e.g., the super- 
symmetric quantum mechanics (SUSY QM) and other equivalent procedures as 
the Darboux transformation, the factorization method, and the intertwining tech- 
Qh| nique [21II3HE]- These constructions have been oftenly applied to Hamiltonians 

having discrete energy levels 1011111313110111111121113111111311121^111111131201 However, 
there are few works involving more general boundary conditions, e.g., on periodic 
potentials for which the physical solutions of the Schrodinger equation have to be 
bounded, and hence the spectrum of the Hamiltonian is composed of allowed energy 
bands separated by the spectral gaps (see e.g. ;21J). Our goal in this paper is to 
show that the SUSY techniques applied to periodic Hamiltonians can produce new 
solvable potentials 1221121111211123112111123 • We will see that when the transformation 
functions are taken as unphysical eigenfunctions of the initial Hamiltonian with fac- 
torization energies inside of the forbidden bands, new periodic or asymptotically 
periodic potentials can be generated [2111201 • The SUSY periodic partners will be 
produced when Bloch transformation functions u(x) are used, while the non periodic 
case will arise for the general u(x) (not necessarily in Bloch form). As a byproduct, 
we will identify an interesting set of Darboux invariant potentials, i.e., those first 
order SUSY partners which become just a displaced copy of the initial potential 

The organization of the paper is as follows. In the second section a brief overview 
of the first and second order SUSY QM will be presented. Then, some simple facts 
about systems involving periodic potential will be discussed. The SUSY QM will 
be then studied as a tool to generate solvable potentials (periodic or asymptotically 
periodic) from an initial periodic Hamiltonian. The mechanism will be applied to 
the Lame potentials, and the paper will end up with a discussion about the Darboux 
invariance. 
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2 Intertwining technique 



Let us consider the following relationship 



HB ] = B j H, (1) 
H = ~ + V(x), (2) 



dx 2 
d* 
dx 2 



H = ~ + V(x) t (3) 



in which the operator intertwines the two Schrodinger operators H, H. Hence, if 
ip n is an eigenfunctions of H with eigenvalue E n , Hip n = E n i/j n , then ip n oc B^ip n ^ 
is an eigenfunction of H with the same eigenvalue, Hip n = E n ijj n . 

In case that £?' is a first order differential operator 

B^ = -^- + a(x,e), (4) 

the Eqs.(|T]|nj) lead to the standard interrelations between a(x, e), V(x), V(x): 

a'(x, e) + a 2 (x, e) = V(x) - e, (5) 
V(x) = V(x) - 2a'(x,e). (6) 

The Darboux formulae are obtained by taking a = [lnu(x)]': 

—u"(x) + V(x)u(x) = eu(x), (7) 
V{x) = V(x) -2[\nu{x)}". (8) 

Thus, a new solvable potential V(x) can be efficiently generated from V(x) if 
one is able to solve explicitly (0) or ((7j) for a certain e, where e is called factorization 
constant because H and H admit the following factorizations: 



H = BB ] + e, (9) 
H = B ] B + e. (10) 

Notice that u(x) have to be nodeless inside the domain of x in order to avoid the 
creation of extra singularities of V(x) with respect to V(x). This immediately leads 
to the typical restriction in the first order case, e < E , where E is the ground state 
energy of H . 

The eigenvalues E n of H for which B<ij) n ^ belong to the spectrum of H. 
The rest of Sp(if) depends on the kernel of B, Bip e = 0. The solution of this last 
equation, ip e oc l/u(x), is an eigenfunction of H with eigenvalue e. According to the 
normalizability of the non-singular ip e , we observe three different cases: 
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If e = Eq and u(x) = if>o(x), it turns out that ip e is non-normalizable =>- 
Sp(if) = {E n , n = 1, 2, . . .} (the level E is 'deleted' in order to get V(x)). 

If e < Eq and non-normalizable f/> £ can be found =^> Sp(i/) = {E n , n = 0, 1, . . .} 
(the strictly isospectral case). 

If e < Eq and normalizable ^ e can be found =>- Sp(iJ) = {e, E n , n = 0, 1, . . .} 
(a level is 'created' at e in order to get V(x)). 



Suppose now that is a second order differential operator: 

By using again (JUIHI), a pair of equations generalizing ((HHEI) are found piTUj: 

"^ + T + i^J +-2~=V(x), (12) 
V r (a:) = y(x) + 2^'(z). (13) 

The solutions of the non-linear second order differential equation f!12l) can be 
found either in terms of those of the Riccati equation (j^J), a(x, ei), a(x, 62), or of 
those of the Schrodinger equation (J7J), U\{x), u 2 (x), E\ 7^ e 2 : 

= 7 ^ ~ ^ V = -pa «*)]'■ (14) 

Contrasting with the first order SUSY, in which «i and u 2 should be nodeless, 
in the second order case the Wronskian W(ui,u 2 ) has to be free of zeros, although 
u\ and u 2 could have nodes. Once again, Sp(if ) depends on the normalizability of 
the two eigenfunctions of H with eigenvalues e\, e 2 which belong to the Kernel of 
B. Their explicit expressions are: 

ip ei oc — -, ip e2 oc — -. 15 

W{ui,u 2 ) W{u 1 ,u 2 ) 

Different cases can be reported: 

• If €1 = Ei, u\ = ipi, £2 = E i+ i, u 2 = ipi+i, the two ip ei , ip e2 are non-normalizable 

Sp(H) = {E , . . . , Ei^i, E i+2 , . . .}. 

• If e\ < e 2 < E and one can find normalizable if) ei , ip €2 => Sp(iJ) = {ei, e 2 , E n , 
" 0.1.2....), 

• If Ei < ei < e 2 < E i+ i and one can find normalizable i/j ei , i/; e2 =>- Sp(if) = 
{Eq, . . . , Ei, ei, e 2 , E i+ i, . . .}. 
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The intertwining technique and the supersymmetric quantum mechanics are 
closely related, namely, the standard SUSY algebra 

[Q,,# ss ] = 0, {Q l ,Q j } = 6 lJ H ss , i, j = 1,2, (16) 

is simply realized by identifying 



Qi = ^(Q f + Q), Q 2 = -) 5 (Q t -Q), 




The first order SUSY QM arises if is the first order operator of (JH), and thus if f 
is linear in the matrix operator H p = 

H ss = (H p -e). (17) 

On the other hand, the second order SUSY (SUSUSY) QM arises if is the second 
order operator (fTTj) . In this case H ss becomes quadratic in H p [HUTU] : 

H ss = (H p -e 1 )(H p -e 2 ). (18) 

In general, in the higher order SUSY QM is an n-th order differential operator, 
n > 1; in this paper we consider just the cases of first and second order. 
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3 Schrodinger equation with 
periodic potentials 

For periodic potentials V(x + T) = V(x) it is convenient to work with the stationary 
Schrodinger equation in the matrix form: 



d 


- 4 - 




dx 


$ _ 





1 

V{x) - E 



(19) 



There is a linear mapping 'propagating' the solution from a fixed point (let us say 
x — 0) to an arbitrary point x: 



i/j(x) 
ip\x) 



b(x) 



f(0) 



(20) 



The 2x2 symplectic matrix b(x) is called transfer matrix. The general behaviour 
of i[) and the spectrum of H depend on the eigenvalues of the Floquet matrix b(T), 
which in turn are determined by the discriminant of b(T), D = D(E) = Tr[b(T)}: 



(3 2 - DP + 1 = 



P± = D/2 ± y/D 2 /4 - 1, (3+(3- = 1. 



(21) 
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The Block functions are particular solutions of ()19|) arising if [ip(0),ip'(0)] T is one of 
the eigenvectors of b(T), i.e., 

^(T)=/3V(0), V'CO=/V(0). (22) 
According to the values of D(E), three different physical behaviours are observed: 

• \D(E)\ < 2. It turns out that j3± G C, |/3±| = 1. This implies that any solution 
(l2llj) is bounded, and hence E belongs to an allowed energy band. 

• \D(E)\ = 2. Both (3+ and /3_ become either +1 or —1, and the associated 
Bloch functions are periodic or antiperiodic respectively. The Floquet matrix 
b(T) is degenerated, and the values of E for which \D(E)\ = 2, denoted as 

E < E x < E v < . . . < Ej < E y < ... 

define the band edges (which belong also to the spectrum of H). 

• \D(E)\ > 2. Now j3± G K, implying that the solutions of (J2*U|) are unbounded. 
Hence, E belongs to a forbidden energy band. 

4 Super symmetrically transformed 
periodic potentials 

It is straightforward to apply the SUSY techniques of section 2 to the periodic 
potentials of section 3 in order to generate solvable potentials from a given initial 
one. Let us employ first as transformation functions the periodic or antiperiodic 
Bloch functions associated to the band edges. The results are the following: 

• 1-SUSY using the 'ground state' eigenfunction ip [SHlElj- It turns out that 
V(x) is non-singular and periodic. As maps bounded eigenf unctions of H 
into bounded ones of H, unbounded into unbounded, etc. H and H have 
the same band structure. 

• SUSUSY employing the band edge eigenfunctions ipj, ipf bounding the energy 
gap (Ej,Eji) The Wronskian becomes nodeless and periodic =^> V(x) is 
periodic. Once again B^ transforms bounded eigenfunctions of H into bounded 
ones of H, etc. =3- H and H have the same band structure. 

A generalization employing Bloch functions with e inside of a forbidden energy 
band can be easily implemented. We distinguish some interesting cases: 

• 1-SUSY using Bloch eigenfunctions for e G (—oo,E ) [22J. Those eigenfunc- 
tions are nodeless, leading to a non-singular periodic superpotential a(x) =^ 
H and H have the same band structure. 
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• SUSUSY employing two nodeless Bloch eigenfunctions associated to 61,62 
which belong to (— 00, Eq) [2E1I2H1- It turns out that (3{x) is non-singular 
and periodic =>- H and H have the same band structure. 

• SUSUSY using two Bloch eigenfunctions with nodes for ei, e 2 such that Ej < 
ei < e 2 < Ej/ Once again /3(x) becomes non-singular and periodic 

H and H have the same band structure. 



a b 




-T T -T T 



c d 




-T t -T T 



Figure 1: The Lame potentials (gray curves) of parameter m = 1/2 and their periodic SUSY 
partners (black curves) generated by employing physical band edge eigenfunctions: a-c) 1-SUSY 
for e = Eq and n — 1,2,3 respectively; d) SUSUSY for e\ — E\, €2 — Ey and n = 3. 

Let us see next what happens when general (non-Bloch) eigenfunctions of H are 
employed as transformation functions [2EH2H]- 

• 1-SUSY using non-Bloch solutions of (j!9|) for e 6 (—00, Eq). It turns out that 
nodeless u can be found such that ip ei is normalizable. The new potential V(x) 
has a local nonperiodicity but it is asymptotically periodic. The operator 
maps bounded eigenfunctions of H into bounded ones of H, unbounded into 
unbounded =>■ H and H have the same band spectrum but there is a bound 
state of H at e. 

• SUSUSY employing two non-Bloch eigenfunctions for 6i,e 2 such that e 1 < 
e 2 < E . ui,u 2 can be selected such that (3 is non-singular and ip ei , ip e2 are 
normalizable. V(x) presents local nonperiodicities but it is asymptotically 
periodic. As maps bounded eigenfunctions of H into bounded ones of H, 
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etc. =>- H and H have the same band structure but there are two bound states 
of H at ei, 62. 

• A similar SUSUSY procedure with e±, e 2 such that Ej < e\ < e 2 < Ey can 
'create' two bound states inside the spectral gap (Ej, Ej'). 

A nice illustration of the SUSY techniques is given by the Lame potentials: 

V(x) = n(n + l)msn 2 (x\m), n e N, (23) 

where sn(x|m) is a standard Jacobi elliptic function of parameter m G [0, 1]. The 
spectrum of H has 2n + 1 band edges defining n + 1 allowed and n + 1 forbidden 
energy bands. Some results applying the SUSY techniques using the band edge 
eigenfunctions are shown in figure 1. Examples employing Bloch functions with 
factorization energies inside of the forbidden energy bands arise in figure 2. Fi- 
nally, some cases implemented by means of non-Bloch transformation functions are 
illustrated in figure 3. 



a b 




-T t -T t 



Figure 2: The Lame potentials (gray curves) with m = 1/2 and their periodic SUSY partners 
(black curves) generated by using unphysical Bloch eigenfunctions: a) 1-SUSY with e = — 1, e < E 
and n = 1; b) 1-SUSY with e = 0.4, e < E and n = 2; c) SUSUSY with e x = 1.6, e 2 = 2.9, 
Ei < ei < e 2 < E v and n = 2; d) SUSUSY with ei = 2.3, e 2 = 5, £1 < d < e 2 < S r and n = 3. 

Let us remark once again that when Bloch transformation functions are used 
(see figures 1-2), the SUSY partner Hamiltonians H and H have exactly the same 
spectrum. However, when non-Bloch solutions are employed (see figure 3), the final 
potential 'acquires' bound states encrusted into the forbidden energy bands, which 
produces local non-periodicities of V(x). These potentials could be useful models 
for the contact effects in solid state physics. 
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a b 




-2T 2T -2T 2T 



c d 




-6T 6T -4T 4T 



Figure 3: The Lame potentials (gray curves) with m = 1/2 and their non-periodic SUSY partners 
(black curves) generated by using unphysical non-Bloch eigenfunctions: a) 1-SUSY with e = 0, e < 
E and n = 1; b) 1-SUSY with e = 0.4, e < £ and n = 2; c) SUSUSY with ei = 1.2, e 2 = 1.3, 
Si < ei < e 2 < £V and n = 1; d) SUSUSY with e x = 1.51, e 2 = 2.51, #i < ei < e 2 < E v and 
n = 2. All the non-periodic potentials have bound states at the positions of the corresponding e's. 



5 Darboux invariant potentials 



Let us notice that for the Lame potential with n — 1 the 1-SUSY technique which 
employs the Bloch functions of either spectral gaps or band edges produces a V(x) 
which is just the initial potential of a displaced argument, V(x) = V{x + 5) (see 
figures la and 2a). That phenomenon was discovered by Dunne and Feinberg for the 
lowest band edge eigenfunction u(x) = ipa(x) with 5 = T/2 (figure la), and those 
potentials were called selfiso spectral (231 ( see a ^ so ESI)- Here we observe a more 
general invariance arising for n = 1 and 5 arbitrary, which is illustrated in figure 
2a. We propose the name translationally invariant under Darboux transformation 
or simply Darboux invariant potentials |2*Hll29| . It would be interesting to seek when 
the 1-SUSY techniques induce that kind of symmetry. The necessary and sufficient 
condition in order that the Bloch transformation function vr{x~) will produce a 
Darboux invariant potential is [231 12*%] 

u p (x)u l/p (x + 5) = constant, (24) 

where w 1///3 (x) is the second Bloch eigenfunctions associated to the factorization 
energy e. The restriction (f2*4*|) is satisfied by the Bloch solutions for the Lame 
potentials (f2*3*|) with n — 1 but it is not for n > 1. We should be able to translate 
the requirement (J24)) into a restriction onto the form of the potential which is going 
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to be Darboux invariant. By using carefully the 1-SUSY techniques it can be shown 
that the Weierstrass potentials are the only Darboux invariant potentials [21]. In 
particular, the Lame potentials with n = 1 are included in the Weierstrass family, 
but there are inside also other interesting non-periodic ones as the 1-soliton well. 
This result explains why the Lame potential with n = 1 is Darboux invariant but 
those with n > 1 are not. It has as well shed some light about the general Darboux 
invariant potentials, not necessarily periodic. 

Acknowledgments. The author acknowledges the support of CONACYT, project 
32086-E. 
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